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ABSTRACT 

Our analysis of the Deep CFHT M33 variability survey database has uncovered 5 Beat Cepheids 
(BCs) that are pulsating in the fundamental and first overtone modes. With only the help of stellar 
pulsation theory and of mass-luminosity (M-L) relations, derived from evolutionary tracks, we can ac- 
curately determine the metallicities Z of these stars. The [O/H] metallicity gradient of —0.16 dex/kpc 
that is inferred from the M33 galacto-centric distances of these Cepheids and from their 'pulsation' 
metallicities is in excellent agreement with the standard spectroscopic metallicity gradients that are 
determined from H II regions, early B supergiant stars and planetary nebulae. Beat Cepheids can 
thus provide an additional, independent probe of galactic metallicity distributions. 
Subject headings: (stars: variables:) Cepheids, stars: oscillations (including pulsations), galaxies: 
individual: M33, galaxies: abundances, galaxies: evolution 



1. INTRODUCTION 

The determination of the abundance distributions pro- 
vides an important constraint for evolution models of 
galaxies. These abundances which are traditionally ob- 
tained through spectroscopic means, although other ap- 
proaches have also been developed, are beset with large 
uncertainties. In this Letter we discuss a totally indepen- 
dent determination of the galactic metallicity (Z) distri- 
bution with the help of Beat Cepheids (BCs hereafter) 
which promises to have a good accuracy. Our approach 
is sensitive to the opacity, and as such it is sensitive to 
the iron group elements, but it is relatively insensitive to 
the chemical makeup of Z , and to Y as well. 

Cepheids are young, intermediate mass, bright, peri- 
odic variable stars. While in the core He burning evolu- 
tionary phase, they lie in an area of the HR diagram, the 
instability strip where their envelopes are vibrationally 
unstable through the kappa mechanism (Cox & Giuli 
1968). They undergo radial pulsations, most of them 
periodic, in a single mode (cither the fundamental mode, 
F, the first overtone, 01, or the second overtone, 02). 
A small fraction of them, called Double-Mode Cepheids 
or Beat Cepheids (BCs), are pulsating simultaneously in 
two radial modes. These pulsations are observed to be 
either a mixture of the fundamental and the first over- 
tone modes F/Ol, or of the first and the second overtones 
01/02. It useful to plot the period ratios vs. period, a 
Petersen diagram (PD) named after (Petersen 1973). It 
is well known that a comparison in a PD of the BCs in 
our Galaxy, LMC and SMC shows that the period ratios 
of F/Ol depend strongly on metallicity (unlike 01/02 
BCs). This behavior has been taken advantage of to con- 
strain stellar models, (e.g., Moskalik, Buchler & Marom 
(1992), Christensen-Dalsgaard & Petersen (1995), Mor- 
gan & Welch (1997)) or to test galactic distances (e.g., 
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Kovacs (2000), Beaulieu, Buchler & Kollath (2001)). 

The M33 late spiral galaxy is known to exhibit a fairly 
steep metallicity distribution. Since Vilchcz et al. (1988) 
various groups have used HII regions to estimate the 
[O/H] distribution in M33. Garnett et al. (1997) found 
an [O/H] gradient (i.e. abundance versus deprojected 
galacto-centric radius) of -0.11 dex/kpc in agreement 
with earlier work. The first studies with B type su- 
pergiant stars yielded a gradient of —0.16 ± 0.06 (Mon- 
teverde et al. 1997), although a more recent estimata- 
tion (Urbaneja et al. 2005) gave —0.06 ± 0.02, which is 
also the value obtained through an analysis of red gi- 
ant branch photometry by Tiede et al. (2004). From 
planetary Nebulae Magrini et al. (2004) derive a metal- 
licity gradient of -0.14. A recent study of HII regions 
by Crockett et al. (2006) however revises the [O/H] gra- 
dient to —0.012 ± 0.011 dex/kpc, much shallower than 
all the other ones. As the data of the above references 
show, one should keep in mind though that representing 
the metallicity distribution as a gradient may be a gross 
oversimplification. 

The M33 variability survey (Hartman et al. 2006) is a 
wide-field (one square degree) survey of this galaxy con- 
ducted with the Canada- France-Hawaii Telescope, in the 
g'r'i' filters, and the first two years of data are available 
with more than 80000 light curves. This survey is par- 
ticularly well suited for variable stars with periods of the 
order of a few days to a few weeks. The excellent rel- 
ative photometry obtained with the image subtraction 
technique (Alard 2000) has given us the opportunity to 
search for BCs at 840 kpc. These BCs allow us then to 
determine the metallicity gradient across the M33 galaxy. 

2. BEAT CEPHEIDS IN M33 

2.1. The Search for Beat Cepheids in M33 

A set of 3023 Cepheid candidates was identified from 
Hartman et al. (2006) and periods were determined si- 
multaneously using the One Way Analysis of Variance, 
AoV (Schwarzenberg-Czerny 1989), as implemented for 
the EROS project (Beaulieu et al. 1997). This phase 
dispersion minimization method is based on strong sta- 
tistical tests and constitutes a powerful way of searching 
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Table 1 
Beat Cepheids in M33 



Table 2 

Metallicities of M 33 Beat Cepheids 



ID 


Number 


«2000 


<52000 


Po 


Pi 


Pio 


ID 


Pio 


dist 


[O/H] 


Z 


[O/H] 


Z 


Z 










(days) 


(days) 








(kpc) (G97) 


(G97) 


(C06) 


(C06) 


from BCs 


A 


121029 


01:34:59.72 


30:52:25.2 


4.70497 3.38510 


0.7195 


A 


0.7195 


3.5 


-0.355 


0.005 


-0.492 


0.003 


0.004 


B 


160520 


01:32:56.82 


30:41:33.8 


3.97755 


2.86107 


0.7190 


B 


0.7190 


3.1 


-0.311 


0.006 


-0.487 


0.003 


0.005 


C 


133449 


01:34:33.43 


30:51:15.6 


3.82707 


2.71407 


0.7091 


C 


0.7091 


2.3 


-0.223 


0.008 


-0.478 


0.004 


0.011 


D 


234922 


01:33:54.63 


30:35:19.8 


6.17640 


4.33313 


0.7015 


D 


0.7015 


1 


-0.080 


0.013 


-0.462 


0.004 


0.0125 


E 


237367 


01:34:03.97 


30:38:08.4 


6.18792 


4.33481 


0.7005 


E 


0.7005 


1 


-0.080 


0.013 


-0.462 


0.004 


0.0135 



Identifications, coordinates, period of fundamental mode (Po); of 
first overtone mode (Pi) and period ratio (Pio). 



Period ratio (Pio); deprojected galacto-centric distance; [O/H] 
and Z metallicities from Garnctt et al. (1997); same for Crockett 
et al. (2006); Z from the BCs. 



for periodic signals of arbitrary shapes, and a confidence 
level is returned for each trial period. Each star has a 
light curve comprising typically 33 measurements in each 
Sloan band g' , r', i'. We chose the r' band data in our 
analysis because of the highest signal to noise ratio. 

For each Cepheid candidate, we fitted a Fourier series 
to the light curve with the most probable period. This 
series was then subtracted and an AoV periodogram of 
this prewhitencd light curve was computed following the 
method described by Bcaulicu et al. (1997). For BCs pul- 
sating in the F and 01 modes the period ratios are known 
to be Pio = Pi/Po ~ 0.7, whereas for those pulsating in 
the 01 and 02 modes it is P\jPi ~ 0.8. We visually in- 
spected the AoV periodograms of both the prewhitened 
and of the full light curves to check for an additional peak 
near the following periods P/0.7, Px0.7, P/0.8, PxO.8. 

In this analysis we uncovered 5 F/Ol type BCs, and 
their properties are listed in Table 1. (We did not detect 
any 01/02 BCs.) The smallncss of this sample reflects 
the difficulty of searching BCs with only 33 measure- 
ments, and the sample is certainly not complete. Ad- 
ditional photometric data (at least 100 measurements) 
would be needed to significantly increase our detection 
rate of M33 BCs. 

We note that these 5 BCs span deprojected galacto- 
centric M33 distances from 1 to 3.5 kpc. In Table 2 we 
also indicate the values of the metallicities they would 
have with the standard choice of the [O/H] metallicity 
gradient -0.11 dex/kpc, e.g., Garnett et al. (1997) and 
with the revised one of -0.012 of Crockett et al. (2006), 
after conversion from [O/H] to [Fc/H]=1.417 [O/H]) as 
in Maciel et al. (2003) and to Z, using Z Q = 0.017. 

2.2. Pulsational Beat Cepheid Models as Function of 
Metallicity 

For the purposes of pulsations a chemically homoge- 
neous stellar envelope is uniquely characterized by three 
astrophysical parameters in addition to its composition 
variables. We will assume here that the helium con- 
tent Y = Y(Z) is a function of Z, as one would expect 
from galactic evolution. We will find that our results are 
relatively insensitive to the exact relation Y(Z). Fur- 
thermore, while the stellar pulsation frequencies are very 
sensitive to the abundance of Fe group metals, they are 
not sensitive that of the remaining elements which we 
assume to be a solar mixture. For computational pur- 
poses the parameters that characterize a Cepheid enve- 
lope are most conveniently taken to be the stellar mass 
M, the luminosity L, the effective temperature T e ff of 



the equilibrium model, and the metallicity Z. The linear 
periods are then functions of these variables only, i.e., 
P k =P k (L,M,T cS ,Z). 

The stellar equilibrium models are constructed and 
their linear (nonadiabatic) periods and growth rates com- 
puted with the Florida pulsation code. Turbulent con- 
vection is approximated with a time-dependent mixing 
length model (Yecko, Kollath & Buchler 1998) and the 
values of the a parameters as defined in Kollath, Buchler, 
Szabo & Szubry (2002) are taken to be {a d = 2.177, a c = 
0.4, a s = 0.433, a n = 0.12, a t = 0.001, a r = 0.4, a p = 0, 
ct\= 1.5}. OPAL (Iglesias & Rogers 1996) and Alexander 
& Ferguson (1994) opacities are used. 

For each BC in our sample we compute a sequence of 
models with given T c s in a reasonable range of values, 
say from 5,300 K to 6,600 K. For each of these T c g we 
iterate on L and M in order to generate the observed 
periods Po and Pi . The acceptable models are those for 
which both the F and 01 modes are linearly unstable. 
We note here that the use of linear periods is adequate 
because nonlinear periods differ at most in the fourth 
decimal figure (e.g., Antoncllo & Aikawa (1998)). 

In Fig. 1 we display for each of the five BCs (sections 
of) the loci of the acceptable stellar models that satisfy 
the observed periods Po and Pi in an M-L diagram. The 
large dots denote the acceptable models that are linearly 
unstable simultaneously in the F and 01 modes. Each 
model along the loci has a different T e jf. The adopted 
composition values, X and Z, are indicated in each sub- 
figure. 

Not all values can represent actual Cepheids because 
the latter move along evolutionary tracks. Stellar evolu- 
tion calculations show that in first approximation these 
tracks are horizontal in an HR diagram, i.e., M = M(L). 
This allows us to narrow down the acceptable BC models. 
In Fig. 1 we therefore superpose a sample of evolution- 
ary M-L relations. The solid curves represent those of 
Bono et al. (2000) as a function of M, Y and Z. The 
short and long dashed curves are derived from Girardi et 
al. (2000) and Alibert et al. (1999), respectively, some of 
them interpolated to the appropriate composition. This 
set of M-L curves should bracket the uncertainties that 
are inherent in the evolutionary track calculations. Note 
that in the M-L diagram, in contrast to the loci the M-L 
curves are relatively insensitive to Z. 

One notes that for each star the width AM (and AL) 
of the models that satisfy the period constraints turn out 
to be narrow and comparable to the uncertainty in the 
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Fig. 1. — Mass-Luminosity plane for 12 different metallicities Z and corresponding X. Each panel shows the possible solutions (thick 
dots) for models with simultaneously unstable F and Ol modes for our 5 BC stars. The curved lines denote various evolutionary M-L 
relations (c/. text).( A color version of this figure is available from the e-paper) 



evolutionary M-L relations. One also notes that as the 
metallicity increases, the curves systematically move to 
the right (to larger M) in a similar fashion for each star. 
The loci first shift very slowly rightward, but once in the 
broad vicinity of the M-L curves, they accelerate and 
cross this region very rapidly, and then slow down again. 
This very fortunate feature allows us to put tight limits 
on the metallicity of each star, despite the spread in the 
M-L relations and despite the spread in T e g. 

From Fig. 1 we can therefore read off the metallicities 
of the 5 stars. Star C, for example, crosses the M-L rela- 
tions when Z=0.011. Considering that for Z=0.010 and 
Z=0.012 star C is comfortably on either side of the pos- 
sible M-L relations, we estimate conservatively that our 



determination of Z is better than ± 0.001. We have ver- 
ified that a small change of Y away from the used values 
has only a minor effect on the loci. The metallicities that 
we infer this way for each star are reported in Table 2. 

Hydrodynamical simulations and the amplitude equa- 
tion formalism demonstrate that beat pulsations can 
occur only in a subregion of the region where the F 
and the 01 modes are linearly unstable (e.g. Kollath, 
Beaulieu, Buchler & Yecko (1998); Kollath & Buchler 
(2001); Kollath, Buchler, Szabo & Szubry (2002)). It 
may therefore be possible to further narrow down the 
metallicity of BCs, but this will require an extended 
survey with nonlinear numerical hydrodynamical simu- 
lations. 



p,/p - ' 

0.74 - 




0.2 0.4 0.6 0.8 

Log P 

Fig. 2. — Petersen diagram for Beat Cepheids pulsating in the 
fundamental and first overtone mode in our Galaxy (open trian- 
gles), LMC (stars), SMC (open circles) and M33 (filled hexagons). 
The linear fits are displayed for our Galaxy (solid), LMC (point) 
and SMC (dash). 



been suggested that stellar rotation may falsify the use of 
the PD as a metallicity indicator as for delta Scuti stars 
(Suarez, Garido & Goupil 2006). We have checked that 
for Cepheid models, and for their lower rotation rates, 
the rotational effect on the PD is negligible. 

Fig. 2 shows that 4 of the 5 M33 BCs have extremely 
large periods when compared to the known MC and 
Galactic BCs. Clearly, the faintness of M33 precludes 
the detection of low luminosity (and thus short period) 
Cepheids. But why are the MCs and the Galaxy devoid 
of such large period BCs? It is not clear at this time 
whether this is a selection effect or an evolutionary ef- 
fect (e.g. the horizontal tracks not extending far enough 
blueward) or a pulsational effect (the occurrence of beat 
behavior is fragile in the sense that it can be easily af- 
fected by small changes in astrophysical parameters such 
as the chemical makeup), or whether some other effect is 
responsible. 



Table 2 also displays the metallicities that are inferred 
for our BCs from the 'standard' metallicity gradient, e.g., 
Garnett et al. (1997). The agreement is found to be 
excellent although star C has a somewhat high Z in our 
study. On the other hand, our results are incompatible 
with those derived from the revised gradient of Crockett 
et al. (2006). 

From our Z values and from the galacto-centric dis- 
tances of the 5 BCS we derive a metallicity gradient 
dLog(Z)/dR = -0.2 dex/kpc which translates to an 
[O/H] gradient of -0.16 dex/kpc that is in good agree- 
ment with the standard determinations, and rules out 
the shallow revised one. 

2.3. Empirical Determination of the Metallicity 

Let us reexamine our derived metallicities in light of 
Fig. 2 in which we have plotted a PD for all the known 
F/Ol Beat Cepheids in the Galaxy, in the LMC (Al- 
cock et al. 1995; Soszyhski et al. 2000) and in the SMC 
(Beaulieu et al. 1997; Udalski et al. 1999). We also show 
the best linear fits for all 3 galaxies (although there is 
no theoretical reason for the constant Z curves to be 
straight lines). Our 5 BCs from M33 are plotted as 
solid hexagons. The 2 stars A and B, with the high- 
est period ratios, lie approximately in the continuation 
of the SMC Cepheids, and indeed their derived metallic- 
ities of Z -0.004 and Z -0.005 agree with that of the 
SMC (Z=0.004). The star C with a derived Z = 0.011, 
falls between the LMC (Z=0.008) and the Galactic BCs 
(Z=0.020). Finally, the remaining two stars, D and E, 
have the lowest period ratios, 0.7015 and 0.7005, and 
with derived Z = 0.0125 and Z = 0.0135, lie between 
the LMC and Galactic BCs. 

It is important to stress that Fig. 2 uses only observed 
quantities, namely periods and is therefore independent 
of any theoretical modeling. With a little act of faith one 
could have 'predicted' the Z of the 5 M33 BCs, although 
for 3 of these 5 stars a huge extrapolation would have 
been required. Our calculations suggest that such an ex- 
trapolation, whether linear of curved, may be warranted. 
However, it will be necessary to make a detailed survey of 
linear and nonlinear BC models to ascertain with what 
precision the mere location of a BC in the PD is able 
to provide an empirical estimate of its metallicity. It has 



To finish it is worth reflecting on why this procedure of 
determining the metallicity of BCs merely on the basis 
of the two observed periods works so well. There are 
two reasons: First, Cepheids evolve along tracks, i.e., 
L = L(M(t),T cS (t), Z(t))\ to the extent that these are 
approximately horizontal we can assume a M-L relation 
in the form L = L(M,Z); second, the subclass of the 
BCs lies in a narrow substrip of the Cepheid instability 
strip {e.g., Kollath et al. 1998, 2001, 2002); thus, for a 
given M and metallicity Z , the effective temperature is 
restricted to lie in a narrow range, T e g « T e g(M, L, Z). 
These two constraining relations, plus those of the two 
measured periods therefore provide the 4 relations that 
are necessary to uniquely specify the Cepheid model (L, 
M, T e ff) and its composition Z . 

3. CONCLUSIONS 

We have searched for Beat Cepheids in the database 
of the Deep CFHT M33 Variability Survey, and we 
have uncovered 5 of them that are pulsating in the 
fundamental and first overtone modes. With the use 
of stellar pulsational modeling and evolutionary M-L 
relations we have derived the metallicities of these BCs. 
From the galacto-centric distance R of the BCs we have 
inferred a 'pulsational' metallicity gradient dLog Z/dr 
= -0.2 dex/kpc that translates to an [O/H] metallicity 
gradient of -0.16 dex/kpc, in good agreement with the 
other determinations in the literature, but it rules out 
the revised value of Crockett et al. (2006). We stress 
that this 'pulsational' method of probing the galactic 
metallicity distribution is independent of spectroscopic 
or other techniques that have previously been used. 

This work is based on observations with Mega 
Primc/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada-France-Hawaii Telescope 
(CFHT) which is operated by the NRC of Canada, the 
Institute National des Sciences de l'Univers of the CNRS 
of France, and the University of Hawaii. We express our 
gratitude to the CFHT staff at Mauna Kea for conduct- 
ing the service mode observations on which this study is 
based. This work has been supported by NSF (AST03- 
07281) at UF. One of us (JRB) has profited from the 2006 
Aspen Center for Physics workshops on galactic evolu- 
tion. He also gratefully acknowledges the hospitality of 
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